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1 Introduction 



CP-violation in the B-system is the subject of an active experimental pro- 
gramme [1-3]. The interpretation of experiments in the framework of the 
Standard Model of particle physics or possible extensions requires substantial 
theoretical input, in particular some knowledge of hadronic matrix elements. 
One important piece in this puzzle is the leptonic decay constant of the B- 
meson, Fb, which in principle can be computed accurately by means of lattice 
QCD. 

Due to recent progress in this field [4,5], a precise computation of decay con- 
stants of light mesons (e.g. Fk) is rather straight forward, at least in the 
quenched approximation [6]. The primary reason is that the renormalization 
of the axial current and the improvement of action and current are known 
non-perturbatively [7-10]. Heavy quarks are more difficult, since their mass is 
large in comparison to the confinement scale and consequently lattices with 
a large number of points are necessary to control lattice spacing effects (see 
e.g. [11,12] for reviews of this problem). As an alternative to the standard 
treatment of quarks as relativistic Dirac fields, approximate treatments, such 
as non-relativistic lattice QCD, have therefore been developed to treat heavy 
quarks. 

On the other hand, in the limit mquark ^ oo, the theoretical treatment sim- 
plifies, as in this limit quarks can be described by the "static effective field 
theory", which is believed to be a local renormalizable quantum field the- 
ory [14,15]. Besides the theoretical interest in this limit of QCD, where new 
symmetries appear [16-18], there are relevant phenomenological applications. 
Computations of quantities such as in the static limit provide important 
checks on the systematic errors made when approximate methods are used. 
Furthermore, a very useful strategy for the computation of B-meson matrix 
elements might be an interpolation between the static limit and quark masses 
of order 1 GeV rather than an extrapolation from light quarks to heavy ones. 
This approach has been tried in the past in order to compute F^ [19,11,12], 
but was limited in precision due to two reasons. First, the renormalization of 
the current in the static effective field theory was known only perturbatively 
and second the computation of the matrix element in that limit turned out 
to be rather difficult. In this paper we assume that the second difficulty can 
be resolved by the methods of [20] or [21] possibly combined with the ideas of 
[22] . For the time being, we concentrate on preparing the ground for removing 
the first obstacle. To this end we introduce a renormalization condition in the 
Schrodinger functional as an intermediate renormalization scheme. The setup 



For a recent review of heavy quark physics on the lattice, we refer the reader to 
[13]. 
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is such that this scheme can be used for a non-perturbative computation of 

the current renormahzation necessary to obtain the renormalization group in- 
variant and scheme independent current, following the strategy described in 
[23]. 

A previous attempt to renormalize the static axial current non-perturbative- 
ly [24] failed when it was applied numerically [25]. We have checked that 
our renormalization condition yields good precision when evaluated in Monte 
Carlo computations [26]. 

Since a systematic discussion of 0(a) improvement of action and axial current 
in the static effective field theory has not been carried out before, we present 
it in Sect. 2. This discussion has some relevance also in the context of non- 
relativistic QCD. 

In Sect. 3 we introduce the Schrodinger functional for static quarks, discussing 
the continuum formulation and the lattice formulation as well as the pertur- 
bative expansion of certain correlation functions. 

In Sect. 4 we then describe the renormalization of the static axial current 
in various renormalization schemes, introducing in particular the SF-scheme 
which is suitable for non-perturbative MC computations. 
The two-loop anomalous dimension of the current in the SF-scheme is pre- 
sented in Sect. 5 and lattice spacing effects in the SF-scheme are computed to 
one-loop order in Sect. 6. After the one-loop computation of c^^* and 6^^*, co- 
efficients necessary for 0(a)-improvement of the axial current, we finish with 
a detailed summary of our results. Appendices contain our notations, some 
details of the perturbative computations and tables with numerical results for 
certain coefficients in the perturbative expansions. 



2 Static quEirks on the lattice 

In this section we discuss QCD with a number A^f of light quarks, which for 
simplicity are taken to be mass-degenerate, and one additional heavy quark 
flavour treated in the static effective field theory . Generalizations are straight 
forward. The theory is formulated on infinite Euclidean space-time, approxi- 
mated by a hypercubic lattice with spacing a. All statements and formulae are 
unchanged if one considers the theory on a torus, imposing suitable periodic 
boundary conditions. 

In order to present a coherent discussion, we start from the well-known ar- 
guments about renormalization of the effective theory and then proceed to 
analyse what is necessary to achieve 0(a) improvement. 
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2. 1 Basic formulation 



The static effective field theory is defined by a path integral over Grassmann- 
valued fields and ip^^ which are associated to the lattice sites. They carry 
colour and Dirac indices. For ease of notation they are taken as 4-component 
spinors, but satisfy the constraint 



P+V'h-V'h, KP+^K. P+ = |(l + 7o), (2.1) 

leaving us with two degrees of freedom per space-time point (the complemen- 
tary components are set to zero). We follow [14] and choose 



= ^'irk^ E V'h(^)(V*o + 5m)M^) , (2.2) 

X 

as our lattice action. The definition of the backward lattice derivative, V*, 
and other conventions are summarized in App. A. The factor 1/(1 -|- a5m) 
is convenient, but irrelevant for physics, since it may be absorbed into the 
normalization of the static fields. The above Lagrangian describes quarks, 
which, being treated in the static approximation, propagate only forward in 
time. 

For completeness we also give the action describing the propagation of static 
anti-quark fields, ■0^, V'ii" 



= «'TTk^ E V'h(^)(Vo + 5m)^l^^{x) (2.3) 

X 

with 



P_V'h = V'h, V'h^'-^V'h, P- = i(l-7o)- (2.4) 

Since the properties of ■0^) '4'h completely analogous to the ones of t/'h, 
we restrict ourselves to the latter from now on. The dynamics of gauge 
fields and relativistic quarks is governed by the action Sq -\- S-p [27,7] listed in 
App. A. 

The time component of the axial current composed of a static quark and a 
relativistic anti-quark field, ■0], assumes a form. 



At\x) = V'.(^)7o75V'h(x) , (2.5) 

familiar from the relativistic case. Observables are defined in the usual way 
through a path integral with total action 
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S — Sf -\- Sq + 5*11 



(2.6) 



2.2 Extra symmetries 



In our discussion of renormalization and 0(a)-improvement we will need two 
symmetries which are present in the effective field theory, but not in finite 
mass QCD. 

1) The static action, eq. (2.2), (as well as the integration measure) is invariant 
under the SU(2) rotations. 



A — VtP^, i^^V-' , (2.7) 

with 

V = exp{-i(f)ieijk(Tjk) , (2.8) 

and transformation parameters (pi. This property is called heavy quark spin 
symmetry [16-18]. 

2) Static quarks propagate only in time. The heavy quark flavour number 
is therefore conserved locally. Mathematically this may be formulated as the 
invariance of action and measure under 



V-h e^''^"^ A, V^h ^he-"^^"\ (2.9) 

with an arbitrary space (but not time) dependent parameter ?7(x).0 The axial 
current carries a flavour number corresponding to the transformation property 



A'^"-' (x) — > e^^'W A'^^' (x) . (2.10) 

As a consequence of the invariance of the action under eq. (2.9), the propagator 
of the static quarks is proportional to a (lattice) 5-function in space (see 
eq. (2.13), below). We may say that the quarks are "exactly static" for flnite 
lattice spacing. 



2.3 Renormalization 



Renormalization of the effective theory has been discussed before [15], also in 
the context of the lattice regularized theory [28-31]. Using the above symme- 
tries, the discussion is quite simple and may be elevated also to the level of 
0(a)-improvement as will be done below. 

^ We thank U. Wolff for pointing this invariance out to us. 
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We point out, however, that one essential assumption has to be made. We 
assume that both Lagrangian and local composite fields can be made finite by 
addition of local counterterms of same or lower canonical dimension. This has 
been shown to be true to all orders of perturbation theory for the relativistic 
theory [32] but not for the effective theory. Of course this is the assumption 
that the static limit is described by a local effective field theory. Perturbative 
computations performed in the past [15,28-30,33] and also the ones presented 
below, provide non-trivial checks that this assumption is justified. 

With this assumption, we note simply that eq. (2.2) contains all local terms 
with mass dimension d < 4 which comply with the standard (lattice-) rota- 
tional symmetry in the three space dimensions as well as the heavy quark spin 
symmetry.[^ Apart from 6m, no counterterm is needed. 

At first sight, the presence of 6m is worrying, as - for dimensional reasons - it 
is linearly divergent, Sm oal/a. Unlike the corresponding mass counter term in 
the relativistic theory, there is no symmetry which allows a non-perturbative 
determination of 5m. Indeed, this represents a considerable problem in the 
determination of the b-quark mass from lattice results in the static approxi- 
mation [34] . It is therefore important to realize that the dependence of correla- 
tion functions on 5m is known explicitly and one can easily find combinations 
where 5m cancels out. 

As an example, we consider the correlation function 

c^AA(xo) = E(^r(^)(^r)^(o)) (2.11) 

X 

with {A^Q^^)'^{x) = '?/'h(x)7o75'?/'i(x). Performing the integration over the quark 
fields (static and relativistic) yields 

Caa{xo) = E( S^{x, 0)S\x, 0))g (2.12) 

X 

with a static propagator. 



Sh{x, y) = 9{xo - yo) 5(x - y) 

X (1 + a5m)-(^«-^°)/'^ W{y, x)-^P+ 

W{x,x) = l (2.13) 

W{x, X + Rfi) = U {x, fi)U{x + afi, fi) x . . . 

xU{x + {R — a)fi, fi) for R>0 



^ The additional flavour symmetry, which is present when several flavours of static 
quarks (e.g. b and c) are considered [16-18], does not play a role in our context. 
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and a rclativistic propagator, satisfying the Dirac equation 

{D+mQ)S{Xj y) = 5{x—y). Here D denotes the Dirac operator, eq. (A. 24). The 
gauge field average (.)g is to be performed with action Sq — log(det(£' + mo)) 
and the trace extends over Dirac and colour indices. The functions 5, 9 are 
defined in eqs.(A.14,A.15). 

If one starts with Sm = (and after renormalizing the mass of the relativistic 
quarks), one discovers that the correlation function eq. (2.11) is linearly diver- 
gent. This divergence is due to the well known self energy of a static quark, 
[29] 



^seif = a-' {eigl + 0{g^)} , ei = ^ x 19.95. (2.14) 

It has to be cancelled by a proper choice of 6m to obtain a finite continuum 
limit of correlation functions in the static effective field theory . The finite part 
of Sm is a matter of convention (scheme) as usual. Just to give an example, 
one could require 

(ao + aS)log(CAA)(a;o) = at = l//x (2.15) 

for some renormalization scale // as a condition to fix Sm. With such a con- 
vention 



(Caa)r(xo) = (Zf^TCAAixo) (2.16) 

is finite for a suitably chosen renormalization constant Zf^^^, since there is no 
other local composite field with dimension d < 3 and the quantum numbers 
of Af'^^ which could induce mixing. 

Alternatively we may use the fact that due to eq. (2.13) 



Caa{xo) = CAA{xo)\sm=o X (1 + a^m)"^"/" (2.17) 

and form ratios of correlation functions where Sm cancels. One then avoids a 
normalization condition such as eq. (2.15). This is quite convenient as also for 
a determination of Fb no knowledge of 6m is required: it may be determined 
from (CAA)R(a;o) "^a;o^oo -^B^^B exp(— £^a;o), where E but not Fb depends on 
Sm. We shall adopt this procedure in our definition of Z^^* in the Schrodinger 
functional in Sect. 3. 
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2-4 0(a) improvement 



We now discuss, which terms have to be added to action and current, in or- 
der to have a lattice formulation where (renormalized) on-shell observables 
approach the continuum with a rate proportional to a^. Such a formulation is 
called 0(a)-improved, as the leading lattice spacing corrections proportional to 
a are removed. We follow the argumentation of [7], where 0(a)-improvement 
is discussed for the theory with relativistic quarks. We assume that the ac- 
tion S'impr, eq. (A. 25), [27,7] is chosen for these fields and consider only the 
modifications necessary in and the static currents. 



2.4.1 Action 

According to Symanzik, the cutoff effects of a lattice theory can be described 
in terms of an effective continuum theory. One starts [7] from the terms which 
contribute to Symanzik's effective action [35]. These may then be cancelled 
by adding terms to the lattice action, resulting in an improved action. 

Symanzik's action contains the lattice spacing as an explicit parameter. It is 
written as 



^eff^S-o + a^i (2.18) 

to the order in the lattice spacing considered. Sq is the standard continuum 

action and Si = J d'^xCi{x), with £1 a linear combination of local fields of 
dimension five with the appropriate symmetries [7] . A basis of such fields is 
given by 



01 = V^h^O^O^h + ^h^O^O^h, 

C»2 = m(^h^o^h - ^h^oV'h), (2.19) 

where denotes the covariant derivative. The mass-parameter, m, is the 
mass of the relativistic quarks. There is no mass in the effective theory (6m 

is a counter term), which could appear here. Furthermore, terms V'h'^oi-^oj'^h 
and ipYi^jkFjkiph do not appear since P^a^jP^ = and since the latter term is 
not spin-symmetry invariant. The reader might also miss a term 

i^^DjDjiP^ + i^^'Dj'DjiPi, (2.20) 
in our list. It is not invariant under the symmetry transformation eq. (2.9). 
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The set of operators, eq. (2.19), can further be reduced [36,7] by using the 
formal equations of motion. 



DoMx)^0, 7jj^{x)Do^O. (2.21) 

They yield immediately that both Oi and O2 vanish on-shell. The field O3 con- 
tributes a mass-dependence to Sm. Such terms play a role in mass-independent 
renormalization schemes [7]. 

In order to cancel the 0(a)-effects, it suffices to modify dm or equivalently to 
add a correction term 



55h = a'J2bh ml M^)M^) , (2-22) 

X 

to the lattice action. The coefficient is a function of the bare coupling qq, 
only; rriq denotes the subtracted bare quark mass [7]. Obviously 6h vanishes 
in the quenched approximation and is of order 0(5^0) in general. 

As explained at the end of Sect. 2.3, for many applications one may consider 
observables where 5m and therefore also foh drop out. 

2.4-2 Axial current 

Similarly to the action, renormalized composite lattice fields are described 
by effective fields in Symanzik's theory. For the time component of the axial 
current we find 

«-*)eff = + a MSA^o^% (2.23) 

k=l 

with some coefficients a;^. A basis for {{SA^^^)k} is 



(Mr)i=?i^,7,75^h, 

(Mr)2 = V'i75/^o^h, (2.24) 

(Mr)3=?i^o75V^h, 

(M'o*''*)4 = "^^i7o75^h. 

Terms containing Dji/j^ do not need to be taken into account since they vi- 
olate eq. (2.10). Requiring improvement to hold only on-shell, this set can 
be reduced: From eq. (2.21) one concludes that {SAq^^)2 vanishes. Due to 
V'h = 7o'0h, the operators {SAf'^^)i, {6Af'^^)s and {SAf'^^)4, are connected by 
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the (Dirac) field equation for the relativistic quark. This means that one of 
these terms can be eliminated. We keep {5Af'^^)i and (5Aq'^*)4. 



In order to achieve a cancellation of the 0(a) lattice spacing effects, we add 
a corresponding combination of correction terms to the axial current in the 
lattice theory. We then write the improved and renormalized current in the 
form 



with a mass-independent renormalization constant Z^'^^ and improvement co- 
efficients, Ca^*, depending again on go but not on rriq. 

2-4-3 Vector current 

The previous discussion carries over to the time component of the vector 
current. Its improved and renormalized version may be parameterized as 



The space components of the vector (axial) current are related by the spin 
symmetry to the time component of the axial (vector) current. Therefore they 
do not need a separate discussion. 

2-4-4 Relation to non-relativistic QCD 

The static effective field theory may be considered to be the infinite mass limit 
of non-relativistic QCD (NRQCD). Also in NRQCD the action is expanded 
in terms of composite fields of increasing dimension. Conventionally the coef- 
ficients are assumed to be proportional to inverse powers of the quark mass 
times "logarithms" . When one wants to formulate non-relativistic QCD includ- 
ing 0(a)-improvement, additional terms proportional to the lattice spacing are 
needed. These are just the terms discussed above. 

The necessity of a term proportional to c^^* has first been noted by Morn- 
ingstar and Shigemitsu [37] in a 1-loop calculation in lattice NRQCD. They 
quote 



(Ag^*)o = (1 + 6f *amq)(Af *)o 
(^f^*)o = + acT'Aljl^i^j + V*)V'i 



'h, 



(2.25) 
(2.26) 




(2.27) 
(2.28) 



(2.29) 
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with 



.stat{l) 



^ X 1.00(1). 



(2.30) 



Our analysis implies that additional terms proportional to the mass of the rel- 
ativistic quarks arc necessary to achieve full 0(a)-improvement when one deals 
with non-zero quark masses in lattice NRQCD. Phrased in the language which 
is usually used, the coefficients in the NRQCD Lagrangian (or fields) have to 
be taken as functions of the NRQCD mass, the bare couphng and aniq. Since 
the dependence on arriq is needed only in the context of 0(a)-improvement, 
the NRQCD coefficients may be Taylor-expanded in that variable. Numeri- 
cally, such terms might be relevant when one deals with the strange quark, 
since NRQCD is usually applied at rather large lattice spacings. 



3 Static quarks in the Schrodinger functional 

In this section, we first introduce the Schrodinger functional in the continuum 
following [38]. As a new element we include static quarks. In Sect. 3.4 we then 
define the Schrodinger functional on the lattice. 

3.1 The Schrodinger functional 

The space-time topology of the Schrodinger functional is a finite cyhnder with 
spatial size L x L x L and with time extent T. In the three space directions, 
periodic boundary conditions are used for the gauge field, whereas the quark 
field boundary conditions are given by eq. (A. 11). Dirichlet boundary condi- 
tions are chosen at the xq — and xq — T boundaries. In the following, all 
space integrals run from to L, all time integrals run from to T. 

For the gauge field Ak{x), boundary fields C(x) and C"(x) are introduced, 
such that 



where P projects onto the gauge invariant content of C"(x) (see [38]). No 
boundary conditions are imposed on the time component of the gauge field. 

Concerning the light quark field, only half of its components are fixed at the 
boundaries. 



Ak{x)\xo=o = C(x), 



A(a:)Uo=T=[lPC'](x), 



(3.1) 
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P+V'i(x) 1^0=0 = pi(x), P_V'i(a:)|^o=r = pi'(x), 

and 



(3.2) 



V^i(x)P_|^o=o = a(x), ^i(x)P+|a;o=r = Pi'(x). (3.3) 

For consistency, the boundary functions must be such that the complementary 
components P-Pi, . . . ,Pi'P- vanish. 

The fermionic action is [39] 

SF[A,i>^,i>i\= j dxo j d^xi>y{x){'yi^Di^ + m)iii{x) 

- J d'^Mx)P+M^)U=T. (3.4) 
3.2 Including static quarks 



As in the case of relativistic quarks, boundary fields ph and p^' defined 
such that 



V'h(a:;)Uo=o = Ph(x), ipM\xo=T^ Ph'i^)- (3.5) 

Note that no projectors are necessary here because of eq. (2.1). For the same 
reason, we have P_^h(a^) = 0. Spatial boundary conditions do not need to be 
discussed since the static quarks do not propagate in space. 

Corresponding to (3.4), a boundary term has to be added to the action: 



5'h[A,V'h,V'h] = y dxo y d^x'ipi^{x)Do'iljh{x) - j d^x [V'h(a;)V'h(a;)]a;o=T- 

(3.6) 

The complete action is then 

5[A V^i, V'l, V^h,^h] = ^g[^] + SF[A,i^„iPi] + 5h[A V^h, V'h], (3.7) 
and the Schrodinger functional is defined as the partition function 
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J D[A] DfV'i] L>[V^i] D[i/Jh] D^^] e-5[^.V'i,V'i,V'h,V'h]^ (3_g) 

which is a functional of the boundary fields. In order to have simple renormal- 
ization properties, expectation values of operators are defined for vanishing 
boundary gauge fields [7] : 

{O) = / D[A] D[^j;,] D[i^,] D[i;^] (3.9) 
Oe-5[A,v;„^.,v;„^H]||_,^^_,^_^,^_^^^^^^^^. (3.10) 

Here, O can — apart from the gauge field and the quark and anti-quark fields 
— contain the light quark "boundary fields" 



and the corresponding heavy ones 



Ch'(x) = — — -, Ch(x) = (3.13) 

()ph (x) dph(x) 



3.3 Some observables 

For later use, we define the correlation functions 

/r (^o) = / dVd^^z {Ar\x)UyhM^)) (3.14) 

and 

/r = / d^ud=^vdVd=^z (Ci'(u)75Ch'(v)Ch(y)75Ci(z)) , (3.15) 
and in particular the ratio 



X{go,L/a)^^S^jEl3l, (3.16) 
which also depends onT/L and 6. 
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Renormalization. 

Wc assume that the Schrodinger functional is finite after adding local coun- 
terterms of dimension o? < 4 to the bulk action, and surface terms of dimension 
d < 3. The first part was discussed already in Sect. 2. In complete analogy 
to the discussion in [40], the surface terms can be shown to be equivalent to 
a multiplicative renormalization of the fermionic boundary fields (i, . . . , Ch', 
viz. 



CiR — -^iCb ■ ■ ■ ) Ch'n — -^hCh' ■ (3-17) 

The mass counterterm {6m) leads to an overall exponential factor in the renor- 
mahzed correlation functions (/a^*)r and (/i*^*)r, which cancels in the ratio 
X. In this ratio, also the wave function renormalization constants cancel, such 
that a renormalized ratio can be written as 



X^ = Zf'-^X. (3.18) 
3.4 Lattice regularization 



The discretization of the Schrodinger functional with nonzero boundary gauge 
fields is described in [38,7]. Here we choose C = C = 0, leading to the 
boundary conditions 



U{x,k)\,,=o = 1, U{x,k)\,^=T = 1 (3.19) 
for the lattice gauge field. The time components remain unconstrained. 

For the quark fields, boundary fields pi,pi,ph and pi',p/,Ph' with the same 
projection properties as in the continuum are defined on the lattice sites with 
Xo = and Xq = T, respectively. The gauge field part of the action, as well as 
the light quark action, is given in App. A. 

Defining 

ip^{x) = if xo < or xo > T, (3.20) 

the lattice action for the heavy quark with Schrodinger functional boundary 
conditions can be written as 

5h[C/,V^h,V'h] = a''Y^Mx)V*oMx)- (3-21) 
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With boundary derivatives defined as in (3.11) and (3.13), we can write down 
a discretized version of our correlation functions, 



/r(^o) = -a'E^«^'(^)Ch(y)75Ci(z)) 



(3.22) 



and 



fi^'--^ E o(^i'N^5Ch'(v)Ch(y)75Ci(z)). 



(3.23) 



u>v,y,z 



Symanzik improvement. 

To cancel 0(a) effects, local operators of dimension 4, which are summed 
over the boundary lattice points at Xq = and Xq = T, have to be added 
both to the light quark action [38,7] and, in principle, to the heavy quark 
action. Writing down all such operators that are invariant under discrete 3- 
dimcnsional Euclidean rotations, it becomes clear that they can be disregarded 
in the static case, as they are zero either because of eq. (2.1) or by formal 
application of the field equations eq. (2.21). 

To achieve complete 0{a) improvement, the axial current improvement term 



(xo) = -a' E o + ^*)^h(a;)Ch(y)75Ci(z)) (3.24) 



has to be added to /a^* (see section 2). We therefore define the (improved) 
ratio. 



Note again that the wave function renormalization constants and the heavy 
quark mass counterterm cancel in this ratio. Apart from its arguments, Xi 
also depends on the ratio T/L and the angle 9. 

3.5 Perturbative expansion 

First of all, we introduce a gauge vector field q^{x), which is the fiuctuation 
around the background gauge field (zero in our case), and write the link vari- 
ables as 



^statystet ^j^jj 




ft\T/2) + a(fl-'ff^\T/2) 



(3.25) 
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U{x,iJ,) = e^{goaqu,{x)}. 



(3.26) 



We then choose the Feynman gauge by the gauge fixing procedure described 
in [8], and add a ghost field term to the action, which however does not 
contribute to our correlation functions at one loop level. 

To achieve a perturbative expansion of X, the correlation functions f^'^^ and 
jstat expanded in powers of the bare coupling, 



/r (^o) = fr^"\xo) + g'ofr^'\xo) + o{g',), (3.27) 

and 



We use the strategy of [8] to compute the expansion coefficients. First the 
Grassmann integrals over the fermionic variables are carried out, and then the 
resulting gauge field dependent observables are expanded in powers oi go- The 
result of the first step is conveniently written in terms of the matrices 



and 



^1 = ^{P+U{X, 0)-'H,{x)},,=T-a: (3.30) 

where i/^i^ci is the classical solution of the Dirac equation, and q is the coefficient 
of the boundary 0{a) improvement term in the action [7]. Corresponding 
matrices are defined for the static quark. 



and 



= ^ Y.{P+U{x, Q)-^H^{x)},,=T-a. (3.32) 
Then the correlation functions can be written as 
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ft\^o) = -l{tT{H,{xy^oH^{x)})G (3.33) 

and 



/f^* = l(tr{ir/i^h})G. (3.34) 

In these formulae, the trace is taken over the Dirac and colour indices, and 
the symbol {■ ■ ■)g means the expectation value with the effective gauge field 
measure (including the fermionic determinant). 

The matrices Hi and Ki are expanded as in [8,41]. The heavy quark matrices 
are written as 



H^x) = Hi'\x) + goHi'\x) + glH^^\x) + 0{gl) (3.35) 

and 

= + 5o4'^ + 5o^f ^ + 0{gl) (3.36) 

with 



^^^ = ^Eni^r(T-a,x), (3.37) 



X 



= ^ E P^{h'^\^) - aqo{x)Hi'\x)U=T-a, (3.38) 



(3.39) 

The exact form of H^^^ , H^^^ , and H^^^ is given in Appendix B.l. Inserting these 
expansions up to 0((?o) into (3.33) and (3.34), we get the expansion of /|'^'^* 
and /i*^* in terms of one- loop diagrams, which can be found in Appendix B.l. 

Using furthermore 

fff = 5ocr ^'Vir + 0{g',), (3.40) 
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with Ca^*^^^ taken from (2.30), the expansion coefficients of Xi, eq. (3.25), 
defined by 

X, = Xi'^ + g',xi'^ + O(^o') (3-41) 
can be calculated. Their exact form is 



xi'\L/a)^^^^^^^, (3.42) 



Stat (0) 



X^'\L/a) = /X*^*^'^(r/2)/V/r*^°^ - |/f^*^°^(r/2) (/^stat (0)^1 -3/2 jstat(l) 

+ a cj^t«/|t^^t(°) (r/2) /^fff^ . (3.43) 



4 The renormalized current in the SF-scheme 



We are now prepared to discuss the renormalization of the current in different 
schemes. In order to separate clearly the dependence on the regularization 
and the scheme dependence due to different normalization conditions, we first 
renormalize the current by minimal subtraction on the lattice. The result- 
ing current has finite matrix elements and can be related to the one in MS- 
normalization by a finite renormalization. We then introduce the SF-schcme, 
which in contrast to the other two schemes is defined beyond perturbation 
theory. 



4-1 Lattice minimal subtraction and MS -normalization 

The ratio X can be made finite by a multiplicative renormalization of the 
coupling, the light quark mass, and the static-light axial current, 

Qr = Zfr{l + bgamq)go, (4.1) 

rriR = Zjy\ + b„,amq)mq, (4.2) 

= Zf''\l + bt^^am^)X. (4.3) 

The parameter mq denotes the subtracted mass of the light quark, 

mq = mo - mc, (4.4) 
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where mo is the bare hght quark mass. At finite lattice spacing, different def- 
initions of the critical mass mc are possible, differing by O(a^) terms. Our 
choice is to introduce the PCAC mass as defined in [7], and then one can 
calculate the bare quark mass at which the PCAC mass vanishes. At each 
order of perturbation theory, this mass can be extrapolated to the continuum 
limit, and it is this limit that we will use as the critical light quark mass in 
our calculations. Note that X is only multiplied by the axial current renor- 
malization constant, as the wave function renormalization constants and the 
mass dependent wave function improvement factors cancel. 

In perturbation theory, coupling, mass and static-light axial current are loga- 
rithmically divergent when taking the continuum limit. We define the lattice 
minimal subtraction scheme by requiring that the coefficients of the renormal- 
ization constants are polynomials in ln(a/x) (without constant parts) at each 
order of perturbation theory, with some renormalization scale /i. 

We chose tuq = mc, which means that we do not have to take care of the mass 
renormalization and the mass dependent improvement terms. Expanding the 
correlation functions up to Qq, only the tree level value Zg — 1 is needed, and 
we remain with the current renormalization constant, which is 

Zf%t = 1 - lolnMgL + 0{gL), (4.5) 
with the scheme independent one loop anomalous dimension [42,43] 



70 = (4.6) 
The minimally subtracted ratio. 



^I,lat = Zf%,Xi , (4.7) 
is finite, and its coefficients. 



XSL = Xi'^ , XSL = - 70 Hat^)Xr' (4.8) 

can be extrapolated to the continuum limit a/L ^ 0. Setting fi = 1/L, and 
using the extrapolation procedure described in [44] we obtained the values 
fisted in Table 1. 

The lattice minimal subtraction scheme is related to other schemes by finite 
renormalization. The most widely used scheme is the MS scheme of dimen- 
sional regularization. In the present context this means that the renormalized 
coupling is defined by the MS prescription, but the renormalization of the 
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)^{0) ^(1) «tat(l) 

^I,lat ^I,lat A A 



0.0 


-1.7320508076 


-0.110337(3) 


0.056(1) 


0.5 


-1.6037987399 


-0.135646(2) 


0.036(1) 


1.0 


-1.4193852065 


-0.147875(1) 


0.016(1) 



Table 1 

Perturbative coefficients for T = L and various 6 in the continuum limit a/L = 0. 

current in the effective theory is fixed by matching to the current in the full 
theory (the normalization of the latter is fixed uniquely by current algebra). 
"Matching" means that matrix elements of the current in the effective theory 
are equal to those in the full theory up to terms suppressed by inverse powers 
of the heavy quark mass. For more details we refer the reader to [45,15]. We 
denote the current normahzed in this way by (^^)o- Its relation to the lattice 
minimally subtracted current, 

(^£f )o = ^l^^r , (4.9) 

is 



(^li)o = (xr)MS,lat«f)0, (4.10) 
(XA^*)MS,lat = 1 + (XA^*)M^,lat^lat + ^(^fat)" 

Here and below, all scale-dependent renormalizcd quantities are taken at the 
same renormalization scale //.. Borrelli and Pittori have computed (Xa^*)SIi + 
by considering matrix elements between quark states, and regulating infrared 
divergences by a gluon mass [33] . Extracting from their results the terms which 
contribute to the (local) current, eq. (2.5), one finds 



(xr)S,at = ^-0-14- (4.11) 

In a separate publication we will report on a computation of the matching be- 
tween effective theory and full theory, in lattice regularization and for infrared 
finite correlation functions (defined in the Schrodinger functional ) [46]. Here 
we already quote a result which is directly relevant in the present context. 



(xr)mM = ^ - 0-137(1) ■ (4.12) 
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It is in agreement with cq. (4.11). The lattice minimal subtraction scheme 
does of course depend on the details of the chosen discretization. Eq. (4.12) 
holds if the improved Wilson action is used for the light quark. Results for the 
unimproved case can be found in [28,29]. 



4-2 The SF-scheme 

The SF scheme is a renormalization scheme using observables defined with 
Schrodinger Functional boundary conditions, and identifying the renormaliza- 
tion scale with the inverse box size L~^. First of all, a renormalized coupling 

is introduced, the definition of which can be found in [47] and the masses of 
the light quarks are set to zero. We then define the renormalized static axial 
current in the SF scheme, 

(AffUii) = Zf% (1 + bt'am^)(Af^% (4.13) 
by using the ratio eq. (3.25), imposing the renormalization condition 

Zf^igo, L/a)Xj{go, L/a) = ^^(L/a) at mo = (4.14) 
and setting = 1/L. 

Perturbatively, the finite renormalization constant connecting the renormal- 
ized currents in the SF and MS schemes can be calculated, 

«^*)o = xT (^^*)o, (4.15) 

xr = i+5isxr^'^+o(5y. (4.16) 

Combining eqs. (4.8,4.10,4.14) and eq. (4.15), and expanding in the coupling, 
we get 



UA ) UA ^MS,lat ;i^(0) ' \^-^' ) 

where it is understood that the last term is to be evaluated for a/L ^ 0. The 
results are contained in Table 1. 
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5 Anomalous dimension in the SF-scheme 



To match perturbative and non-perturbative results with sufficient precision, 
one needs the static-light axial current's anomalous dimension 7(^), defined 
by the renormalization group equation 



/^^^^ = 7(^)(^r)o, (5.1) 

at two loop order in the SF scheme. We calculate the coefficient 7p in the 
expansion 



7'^(^sf) = -^iF(7o + iffsF + ■■■) (5.2) 
by conversion from the MS scheme to the SF scheme. 

The two loop anomalous dimension in the MS scheme is known from [45,48,49], 



1 1 97 

Here C is the Riemann zeta function, and A^f is the number of relativistic quark 
flavours. 

The renormalized couplings in the two schemes are related by 



were the 1-loop coefficient is 



^g 



^^^/, = -^(c,o + c,,iVO (5.5) 



with [8,50] 



ci,o = 1.25563(4), ci,i = 0.039863 (5.6) 

for the SF-coupling defined at ^ = 7r/5. The relation between the currents in 
the two schemes is known from Table 1. 

Analogous to [51], where the discussion is carried through for the anomalous 
dimension of quark masses, we have 
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7r=7r+2&oxr^^-7ox(^\ (5.7) 

where 

is the universal one loop coefficient of the (5 function. 

From our result for (xa'^*)^^'* the T — L case, Table 1, we get 



jf{e = 0.0) = ^{0.52(2) - 0.0733(13)7Vf}, (5.9) 

^f{9 = 0.5) = .^{0.08(2) - 0.0466(13)7Vf}, (5.10) 

jf{e = 1.0) = -^^{-0.36(2) - 0.0199(13)Arf} . (5.11) 

For the range of 9 considered, these coefficients are rather small and the per- 
turbative series appears to be well behaved. 



6 Lattice spacing effects in the SF-scheme 



The renormalization constants at energy scales 1/L and 1/{2L) are connected 
by the step scaling function, 



This lattice step scaling function has a continuum limit (7^^*, 



Ef * {u, a/L) = aT' {u) + 0(a/L) . (6.2) 

With Ua'^*, a non-perturbative renormalization group can be constructed with- 
out taking reference to the lattice regularization. 

As Sj^*(-u,a/L) will be the central object of non-perturbativc studies, it is 
interesting to investigate its discretization errors at the one loop level. As an 
estimate for these errors, we have calculated 



S(u, a/L) . "^'"'"'^itl;/"'"' = + + . . . (6.3) 
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Fig. 1. Discretization errors in the 0(a) improved step scaling function at one loop 
order of perturbation theory. Points are for L/a = 4, 6, 8, 

for three values of 9. The results are displayed in Fig. 1, showing that lattice 
spacing effects - which are 0((a/L)^) since we use 0(a) improvement - are 
moderate at one loop order. Again, we have used rric extrapolated to a — 
in this exercise. In MC computations one has to use a definition of rric at 
finite a/L instead. Then lattice spacing errors are somewhat different. We 
have checked that in the case at hand these differences are tiny when typical 
choices for rric, such as the one in [52], are made. 



7 The improved eixial current at 1-loop order 



Throughout our calculations, wc have used the value from eq. (2.30) for the 
current improvement coefficient at one loop order, obtained by Morningstar 
and Shigemitsu [37] using NRQCD methods. We have also computed this 
coefficient from our data for 6 = 0.5 and 6 = 1.0. First of all, we calculated 
the ratio X without the current improvement term, and then extracted its 
0(a) part by the fitting procedure described in [44], using the formula given 
in Appendix B.2. The results are 



47r 



X 1.01(7) at e = 0.5 



(7.1) 



and 



^stat{l) 



-— X 1.01(5) 



at 0^ 1.0, 



(7.2) 
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in agreement with eq. (2.30). We note that our values also agree with the more 
recent and precise calculation of [53], 

^jat{i) ^ _ J_ ^ 1.0351(1) 

From a similar analysis of data obtained with a nonzero light quark mass, the 
coefficient 6^^' introduced in Sect. 2.4 can be calculated at one-loop order. 
Our result, computed from eq. (B.19) for ^ = and z = Lm\sx = 1, is 

hf^' = I + ht'^'^^^gl + . . . , hT^'^ = -0.056(7). (7.4) 

Here the error is due to the uncertainty in the a/L ^ limit. Determinations 
with other values of z and 9 turned out to be entirely consistent with eq. (7.4). 



8 Summciry 

Our paper covers two main topics. The first one concerns the formulation of 
the static quark effective theory on the lattice, in particular the question of 
0(a) improvement. The second topic is of a more practical nature and deals 
with the formulation of a renormalization condition for the static light axial 
current which can be used to renormalize the current non-perturbatively by 
MC computations. 

I. We have started from the general assumption that the divergence structure 
of static quark effective field theory can be inferred from simple dimensional 
counting and that the theory can be made finite by addition of local countert- 
erms. 

Clearly it would desirable to have a proof that this assumption is correct to all 
orders of perturbation theory. In absence of a proof, we may, however, perform 
checks through perturbative computations as well as finally MC computations. 
All results of our one loop computations are in agreement with the assump- 
tion. Particularly convincing is the agreement of the renormalization factor 
eq. (4.11) of [33] with our result eq. (4.12), obtained from a completely dif- 
ferent quantity. A similar statement holds for the renormalization beyond the 
leading order, i.e. 0(a) -improvement: the coefficient CyJ^* from [37,53] yields 
0(a)-improvement for our correlation functions. 

In more detail, the following terms are found to be necessary and sufficient 
for renormalization and improvement of the theory. 
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(7.3) 



o In the Lagrangian only a mass counterterm dmipY^{x)tpi^{x) is necessary at 
the leading order in the lattice spacing (i.e. to obtain a finite continuum 
limit). 

o The static light axial current is renormalized multiplicatively. 

o For vanishing light quark masses, the static sector defined with the action 
of Eichten and Hill, eq. (2.2), is 0(a)-improved. At finite quark mass, a 
correction term SS^^ = a^J2x bh'>Ti^ip^{x)iphix) has to be added to achieve 
0(a)-improvement. Just like Sm, such a term affects correlation functions 
only in a trivial manner. 

o Improvement of axial and vector current necessitates two 
improvement terms each. For example the renormalized axial current may 
be written as 

(^r )o = (1 + bTam^) [AT' + acrfe75|(V, + V*)t^h 

with Z^^^ the renormalization constant defined at vanishing quark mass and 
two improvement constants Ca^*,&a^*. In perturbation theory, we have 

6ft = l + 6ft« + 0((7o^) 

^jat^^stat(l)^2^0(^o')- (8-1) 

Indeed, Ca^*^^-* had already been determined in the hterature [37,53] and 
was confirmed once again by us. It is rather small. The one loop coefficient 

^stat(i) jg gjyg^ (7.4). 

o Conveniently, the renormalization of the Schrodinger functional with static 
quarks needs only a multiplicative renormalization of the "boundary 
fields" . No new 0(a) (surface-) counterterms appear when the light quark 
mass is set to zero. 



Apart from the special situation of the Schrodinger functional , the renormal- 
ization properties of the effective theory have been known already and are 
included above only for completeness. On the other hand, the general struc- 
ture of 0(o)-iniprovcmcnt has not been derived before. We remark that terms 
such as b%^^ are necessary also in NRQCD, see end of Sect. 2. 

II. We have introduced a renormalization condition for Z^^^ in the framework 
of the Schrodinger functional . This is done such that a non-pcrturbative 
computation of the scale dependence as well as the value of Z^'^^ at a particular 
value of the renormalization scale are possible, in complete analogy to the case 
of the QCD quark masses [52] . 

o In order to be able to follow the strategy of [52], the two loop anomalous 
dimension has been computed; it is given at the end of Sect. 4. 

o Lattice spacing effects in the step scaling function (describing the scale de- 
pendence of the current) have been computed to one loop. They are small at 
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this order, promising the possibihty of controlled continuum extrapolations 
of future MC results. 

MC computations of the step scaling function are in progress [26]. 

As a future perspective, one can follow the same kind of strategy also to renor- 
malize the operators responsible for BB mixing in the static approximation. 
While this is very attractive in particular in view of the phenomenological 
relevance, one should be aware that new difficulties will appear: in this case 
there is a scale-dependent mixing and more terms will be necessary to achieve 
O (a)-improvement. 

Acknowledgement. We would like to thank Peter Weisz, UUi Wolff and 
Stefan Sint for useful discussions. Particular thanks go to Martin Liischer 
who's FORTRAN programs [8] were the basis for our programs computing 
the one-loop correlation functions. 



A Notations and action for relativistic queirks and gauge fields 

We use the conventions of [7] . Together with a few extensions they are hsted 
here for the readers convenience. 



A.l Index conventions 

Lorentz indices /i,!/, . . . are taken from the middle of the Greek alphabet and 
run from to 3. Latin indices k,l, . . . run from 1 to 3 and are used to label the 
components of spatial vectors. For the Dirac indices capital letters A,B, . . . 
from the beginning of the alphabet are taken. They run from 1 to 4. Colour 
vectors in the fundamental representation of SU(A^) carry indices a, (3,... 
ranging from 1 to N, while for vectors in the adjoint representation, Latin 
indices a,b, . . . running from 1 to N'^ — 1 are employed. 

Repeated indices are always summed over unless otherwise stated and scalar 
products are taken with euclidean metric. 



A. 2 Dirac matrices 

We choose a chiral representation for the Dirac matrices, where 
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The 2x2 matrices are taken to be 



Co = -1, Ck = -iak, (A.2) 
with (Tfc the Pauh matrices. It is then easy to check that 

ll = lix, {lix, lu} = 25^^. (A.3) 
Furthermore, if we define 75 = 7o7i7273, we have 



75= (i °). (A.4) 



,0 -1 

In particular, 75 = 7I and 75 = 1. The hermitean matrices 



are exphcitly given by 



-a,;' -^no ak 

where eijk is the totally anti-symmetric tensor with £123 = 1 
A.3 Lattice conventions 



Ordinary forward and backward lattice derivatives act on colour singlet func- 
tions f{x) and are defined through 



dj{x)^-[f{^+af^)-mi 

a 

dlf{x)^^[f{x)-f{x-ail)], (A.7) 

where p, denotes the unit vector in direction /i. The gauge covariant derivative 
operators, acting on a quark field ip{x), are given by 
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V^ipix) = -[X^U{x,n)ij{x + afi) - i>{x)], 

W%ip{x)^-[il){x) - X~^U{x - afi, ii)~^ip{x - aft)], 
a 

The phase factors 



(A.8) 
(A.9) 



y^^^^iae./L^ ^0 = 0, -7r<^fe<7r, (A. 10) 

are equivalent to imposing the generahzed periodic boundary conditions 



il^ix + Lk) = g'^'"iI){x), ilj{x + Lk) = ilj{x)e' 



-idk 



(A.ll) 



They depend on the spatial extent L of the lattice and are all equal to 1 on 
the infinite lattice. The left action of the lattice derivative operators is defined 

by 



^p{x)V 



1 r 

a 
1 



il){x)V^ = - il){x) - il){x - ajl)U{x - ajl,ij)\^ 



Our lattice version of (5-functions are 



(A.12) 
(A.13) 



3 3 

5{x^) = a-^5^^Q , 5(x) = JJ 5{xk) , S{x) = S{x^) 

k=l /i=0 

and we use 



(A.14) 



9{Xf,) = 1 for > 
0{xij) — otherwise 



(A.15) 



A. 4 Continuum gauge fields 



An SU(A^) gauge potential in the continuum theory is a vector field ^^^(a;) 
with values in the Lie algebra su(A^). It may thus be written as 



(A.16) 



with real components Al^{x) and 
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(A.17) 



The associated field tensor, 



(A. 18) 



may be decomposed similarly and the right and left action of the covariant 
derivative is defined by 



The abelian gauge field i9^/L appearing here corresponds to the phase factors 
A^, eq. (A. 10), in the lattice theory. 

A. 5 Lattice action 

Let us first assume that the theory is defined on an infinite lattice. A gauge 
field U on the lattice is an assignment of a matrix U{x,n) G SU(iV) to every 
lattice point x and direction n — 0,1, 2, 3. Quark and anti-quark fields, il^{x) 
and ip{x), reside on the lattice sites and carry Dirac, colour and flavour indices. 
The (unimproved) lattice action is of the form 



where Sq denotes the usual Wilson plaquette action and Sp the Wilson quark 
action. Explicitly we have 



with go being the bare gauge coupling and U (p) the parallel transporter around 
the plaquette p. The sum runs over all oriented plaquettes p on the lattice. 
The quark action. 



D^ij{x) = {df, + A^ + ie^/L)^{x), 
i^{x)D^^i^{x){d^-A^-i9jL). 



(A.19) 
(A.20) 



(A.21) 



(A.22) 



S^[U, ^, ^ ] = ^ ij{x){D + mo)^(x). 



(A.23) 



X 



is deflned in terms of the Wilson-Dirac operator 



\ {7.(V* + V,) - aV* V,} , 



(A.24) 



30 



which involves the gauge covariant lattice derivatives and V*, eq. (A. 7), 
and the bare quark mass, mo- 

The improved action is given by [27,7] 



SunAU, ^] = S[U, ^, ^] + 5S\U, ^, V], (A.25) 
5S[U,i^,ij]=a''Y.c,^i^{x)\a^,F^,{x)ij{x), (A.26) 

X 

with 

F^Ax) = ^ {Q^,{x) - Q,^{x)} , (A.27) 

Qnv{x) — U {x, n)U{x + afi, v)U{x + az>, fi)~^U{x, u)'^ (A. 28) 

+ U {x, v)U{x — ajl + au, n)~^U {x — afi, u)~^U {x — aft, ji) 
-\-U{x — afi, IJ>)~^U {x — afi — au, 

xU {x ~ afi — aO, ^)U{x — au, v) 
+ U{x — av, i')~^U {x — av, ii)U{x + afi — au, i')U{x, 

Turning now to the case of Schrodinger functional boundary conditions, the 
pure gauge action on the lattice is [38] 



Sg[U] = \j2Hp)tr{l - U{p)}, (A.29) 
9o p 

with a weight w{p) equal to 1 for all p except for the spatial plaquettes at 
Xq — and Xq — T, where w{p) — |. 

The quark action in the Schrodinger functional can be written exactly as in 
eq. (A. 23) when the conventions 



i/j{x) = if xo < or xo > T (A.30) 

and 

P-ipix)\,,=o ^ P+ip{x)Uo=T ^ (A.31) 
are adopted. 



B Static quEirks in perturbation theory 
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B.l Heavy quark Feynman rules 



Using the lattice action (3.21) for the heavy quark, we get the heavy quark 
propagator 



S^ix, y) = U{x - aO, 0)"^ . . . U{y, 0)-'e{xo - 2/o)5(x - y)P+, (B.l) 
with 9{xo) and 5(x) defined in (A.15,A.14). Its fourier transform is defined 

by 



S^i^^y) = jsJ2^'''^'"''^S^i^o,yo;p). (B.2) 

Expanding in terms of go, we find 



S^^\xo,yo;p)^P+. (B.3) 
The classical equation of motion for the heavy quark, 

VUhA^) = 0, (B.4) 
is solved by 

V'h,ci(a^) = U{x - aO, 0)-^ ...U{x- XqO, O)-Vh(x), (B.5) 
which leads to 

Hh{x) = U(x - aO, 0)-^ ...U(x- XqO, 0)-^P+. (B.6) 
Expanding this in powers of the bare coupling, and writing 

?o(x) = L-3^e^P-?~o(xo,p) (B.7) 
p 

we get 

Hi'\x) = P+, (B.8) 
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Xo=0 
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Xo=0 



T/2 



2a 2b 3 

Fig. B.l. One loop diagrams contributing to f^^^{T/2). 



T-a 



(B.9) 



p uo=a 



and 



(B.IO) 



with 



(^)i = 7^ E E E e'^^^^^-si'\xo, uo; q+ p)goK - a; q) (B.ll) 

q,p uo=a vo=a 



xSl '{uo - a,vo;p)qo{vo - a;p)P^ 



(B.12) 



and 



,2 T-a 



hI'\x), = — Y: E e^('»+P)-^f (a;o,«o;0)go(«o-a;q) 



ore 

^-^ q,p wo=a 



(B.13) 
(B.14) 



Contracting the gauge field fiuctuations q then gives the diagrams in Fig. B.l 
and Fig. B.2. 



B.2 Calculation of {(f^^^Y^^ and {ff^^^Y^^ 



For the calculation of {c^^^^)^^\ the one-loop coefficient of the ratio Xi is de- 
composed as 
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Xo=0 



.o. 

T-a T 



6b 






Xo=0 



T-a T 




Fig. B.2. One loop diagrams contributing to ff'^'^. The dotted lines denote the link 
from T — a to T. 

xi'\L/a) = X«(L/a)|e,=i + cf ^a4°^(L/a) + ct'^'^aX^'^- (B.15) 
Defining the derivatives d and d* by 



9/(L/a) = l{/(^)-/(L/a)} 
d*f{L/a) = Uf{L/a)-f{i^)}, 

(Jj 

the 0(a) part oi X^^ can be extracted to get 



(B.16) 
(B.17) 



,stat(i) lim„/L-.o g(a + y ) |c-.=i - lim„/^->o c^^ LXi'\^J 

= . „rm.r. l^-^^J 



with the hght quark mass set to zero. 

In a similar way, (6^^*)^-'^) can be calculated, 



istat(l) 



I lim f (a + a*)xW(^)|a,=i - lim cf^LX^'V^^ 
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where 



XW = (1 + 6f t(°)am(°))X(i) + (Zl^Ji'^ + 6f *^°^am«)XW. (B.20) 
The quantity 

^ = ■^m,iat(l + &marnq)Lmq (B.21) 

is kept fixed when taking the hmit. The renormahzation and improvement 
constants for the hght quark mass are known from [41]. 

C Numerical results 



The following results were obtained at T = L, with the bare light quark mass 
set to the critical mass 



with [8] 

amW = -0.2025565 x Cp- (C.2) 
The improvement coefficients used are 

Csw — 1) (C-3) 

and 

Q = l + cJSo' + 0(^o'), (C.4) 
with [41] 

cP = -0.01346 X Cp. (C.5) 

The tree level results can be obtained analytically. Using the notation of [8], 
they are 
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-((M(p+)+Po)e-'"(P^)('^-^°)} at p = 0, (C.6) 



and 



-((M(p+)+Po)e-'"(P^)(^+")} at p = 0. (C.7) 

Also the axial current improvement term can be calculated analytically at tree 
level, 



/|^-t(°)(a;o) = -^5^{e-"(P^)^°-e-'^(P^)(2^-^°)} at p = 0. (C.8) 

K[p ) 

The one loop results are given tables (C.l)- (C.3). They have been obtained 
in double precision calculations. Where an error is given, it has been taken 
from the comparison with the corresponding quadruple precision result. 
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L/a 




jstat(l)(^/2) 




jstat(l) 


4 


1, 


,1696384959451731(7) 


-2, 


,788717306605536(2) 


6 


1, 


,5721721053841535 


-3, 


,744847208299456 


8 


2, 


,049414037209983(4) 


-4, 


,76922384055426(4) 


10 


2. 


,543921176114074 


-5. 


,80311014058304 


12 


3. 


,04401042107373(4) 


-6. 


,8365479914505(2) 


14 


3. 


,54641935327602 


-7. 


,86798061045362 


16 


4. 


,04997105269991(6) 


-8. 


,897379383110(2) 


18 


4. 


,55415464879211 


-9. 


,924989104618 


20 


5, 


,058717578768(1) 


-10, 


,951072803831(9) 


22 


5, 


,5635226305105 


-11, 


,975859714707 


24 


6, 


,0684897524216 


-12, 


,999539542806 


26 


6. 


,5735695645449 


-14. 


,02226731358 


28 


7. 


,0787302119581 


-15. 


,04416980198 


30 


7. 


,583950374904 


-16. 


,06535128339 


32 


8. 


,089215337913 


-17. 


,08589819107 


34 


8. 


,594514673902 


-18. 


,10588276384 


36 


9, 


,099840825103 


-19. 


,12536587521 


38 


9, 


,60518820413 


-20, 


,14439922726 


40 


10. 


,11055260756 


-21, 


,1630270594 


42 


10. 


,61593082380 


-22. 


,1812874905 


44 


11. 


,12132036473 


-23. 


,1992135777 


46 


11. 


,62671927843 


-24. 


,2168341613 


48 


12. 


,1321260162(3) 


-25. 


,234174544(1) 



Table C.l 

f^^^ and /f*^* at one loop order with 6 = 0.0 
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L/a 




jstat(l)(^/2) 




jstat (1) 


4 





.707462157465912(1) 


-1, 


.828520866963525(4) 


6 


1 


.0665692749675324 


-2, 


.53494068630446456 


8 


1 


.45264489740883(2) 


-3, 


.268670455891038(1) 


10 


1, 


.845191095078847 


-4, 


.0063158459818604 


12 


2, 


.23966783029017(2) 


-4, 


.7434246150126(1) 


14 


2, 


.634843388684698 


-5, 


.47929315793510 


16 


3 


.03031175469344(4) 


-6, 


.2139336444356(4) 


18 


3 


.42591623855986 


-6, 


.9474888408264 


20 


3 


.82158857585887 


-7, 


.6801091376222 


22 


4, 


.21729615057110 


-8, 


.4119266083201 


24 


4, 


.61302225619736 


-9, 


.1430518919882 


26 


5, 


.00875788592884 


-9, 


.8735765450344 


28 


5, 


.40449800734282 


-10. 


.603576331834 


30 


5, 


.80023974629234 


-11. 


.333114250324 


32 


6, 


.19598144778345 


-12. 


.062243035152 


34 


6, 


.59172216643778 


-12. 


.791007154176 


36 


6, 


.98746137900958 


-13. 


,51944438137 


38 


7, 


.3831988169726 


-14, 


,24758703382 


40 


7, 


.7789343665427 


-14, 


,97546294777 


42 


8, 


.1746680077261 


-15. 


.70309625284 


44 


8, 


.5703997766495 


-16. 


.43050799027 


46 


8, 


.9661297419670 


-17. 


.1577166104 


48 


9, 


.36185799001(1) 


-17. 


.8847383757(5) 



Table C.2 

f^^^ and /f*^* at one loop order with 6 = 0.5 
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L/a 


j|tat{l)(^/2) 


ystat{l) 


4 


0.3044333805556540(7) 


-0.869671412347571(1) 


6 


0.5486322003500290 


-1.1948835905593809 


8 


0.793695141003886(4) 


-1.53524521245116(1) 


10 


1.038456376564977 


-1.881017287644124 


12 


1.28269652915735(3) 


-2.22894706469196(3) 


14 


1.52651219643810 


-2.57783369422662 


16 


1.77002559881285(11) 


-2.92715777151496(7) 


18 


2.01332870498098 


-3.27666836444461 


20 


2.2564839354402(4) 


-3.6262359702708(2) 


22 


2.4995329275626 


-3.97579139115074 


24 


2.7425037517644 


-4.3252973303331 


26 


2.9854157215170 


-4.6747341210386 


28 


3.2282824349472 


-5.0240921351849 


30 


3.4711136865628 


-5.3733675663630 


32 


3.7139166829809 


-5.7225600072263 


34 


3.9566968307148 


-6.071671021232 


36 


4.199458257584 


-6.420703284271 


38 


4.442204165276 


-6.769660062007 


40 


4.684937072624 


-7.118544889597 


42 


4.927658986791 


-7.467361375608 


44 


5.170371525823 


-7.81611308330 


46 


5.413076007854 


-8.16480346057 


48 


5.65577351702(2) 


-8.51343580067(7) 



Table C.3 

f^^^ and /f*^* at one loop order with = 1.0 
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